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Abstract 

We perfect the recursion-transform method to be a complete theory, which can derive the 
general exact resistance between any two nodes in a resistor network with several arbitrary 
boundaries. As application of the method, we give a profound example to illuminate the usefulness 
on calculating resistance of a nearly mxn resistor network with a null resistor and three arbitrary 
boundaries, which has never been solved before since the Greens function technique and the 
Laplacian matrix approach are invalid in this case. Looking for the exact solutions of resistance 
is important but difficult in the case of the arbitrary boundary since the boundary is a wall or 
trap which affects the behavior of finite network. For the first time, seven general formulae of 
resistance between any two nodes in a nearly mxn resistor network in both finite and infinite 
cases are given by our theory. In particular, we give eight special cases by reducing one of general 
formulae to understand its application and meaning. 
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1. INTRODUCTION 


The computation of the two-point resistance in networks is a classical problem in circuit 
theory and graph theory since the German scientist Kirchhoff found the node current law 
and the circuit voltage law in 1845 [1]. Modeling the resistor network has become a basic 
method to solve a series of complicated problems in the fields of engineering, science and 
technology. Such as the computation of resistances is relevant to a wide range of problems 
ranging from classical transport in disordered media [2], first-passage processes [3], lattice 
Greens functions [4], random walks [5], resistance distance [6,7], to hybrid graph theory[8]. 
For an up-to-date list of references, see Ref. [9]. Modeling the resistor network can also help 
to carry on the research of exact finite-size corrections in critical systems [10-15]. However, 
it is usually very difficult to obtain the exact resistance in a complex resistor networks with 
arbitrary boundary since the boundary is a wall or trap which affects the behavior of finite 
network. It requires not only the circuit theory but also the innovation of mathematical 
theories and methods [16-27]. The construction of and research on the models of resistor 
networks therefore make sense for theories and applications. 

For the explicit computation of two-point resistances in a resistor network, Cserti [16] 
evaluated the two-point resistance using the lattice Greens function. His study is confined 
to regular lattices of infinite size. In 2004 Wu [17] formulated a different approach and 
derived an expression for the two-point resistance in arbitrary finite and infinite lattices in 
terms of the eigenvalues and eigenvectors of the Laplacian matrix. The Laplacian analysis 
has also been extended to impedance networks after a slight modification of the formulation 
of [18]. The Laplacian method has solved a variety of types of resistance network with the 
normal boundary, besides the previous cases [17], there are also the recent Tan-Zhou-Yang 
conjecture [19], and a globe [20]. But the applications of the Laplacian approach require 
a complete knowledge of the eigenvalues and eigenvectors of the Laplacian straightforward 
to obtain for regular lattices. However, the actual computation depends crucially on the 
geometry of the network and, for nonregular lattices such as a network with an arbitrary 
boundary, it is difficult to solve the eigenvalue problem with arbitrary elements. 

In recent years, a new method is created by us [21,22], which is easy to solve the problem 
of nonregular lattices of cobweb and globe[22-25]. This methods splits the derivation into 
three parts. The first creates a recursion relation (expressed by matrix) between the current 
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FIG. 1: An m x n resistor network with an null resistor and three arbitrary boundaries. Bonds 
in the horizontal and vertical directions are resistors r and tq except for four boundaries. 

distributions on three successive vertical lines. The second part diagonalizes the matrix 
relation to produce a recurrence relation involving only variables on the same transverse 
line. The third part derives a recursion relation between the current distributions on the 
boundary. Basically, the method reduces the problem from two dimensions to one dimension. 
Here we refer to this method as the recursion-transform approach ( call it R-T method ). 

Using the R-T method to compute the equivalent resistance relies on just one matrix 
along one direction, and the resistance is expressed by single summation [21-27]. Especially, 
the R-T method is always valid to calculate the resistance when two opposing boundaries 
change. For example, the R-T method is used solving a cobweb network in 2013 [21, 22], 
next a globe network [23], and a fen and a cobweb network [24], very recently a cobweb 
network with 2r boundary is solved [25], and a hammock network with free boundary is 
solved[26], and a resistor network with an arbitrary boundary is solved by the R-T method 
[27] . 

However the original R-T method is still in its infancy, which has to be perfected. In 
real life the boundary conditions of the resistor network are often complicated. This paper 
focus on perfecting the R-T method to be a complete theory, which can derive the gen¬ 
eral resistance between any two nodes in the nearly m x n resistor network with arbitrary 
boundaries. As applications of the R-T method, we gives three examples to illuminate the 
usefulness on calculating resistance of a resistor network, and seven general formulae of the 
resistance between any two nodes in an m x n resistor network in three cases are given. It 
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shows the R-T method is very interesting and useful to solve the problem which have never 
been solved before. 

2.A GENERAL THEORY AND APPROACH 

Fig.l is called a nearly rectangular m x n resistor network with an null resistor and three 
arbitrary boundaries, where r and r 0 are, respectively, resistors in the horizontal and vertical 
directions except for four boundaries, and where m and n are, respectively, the number of 
meshes along vertical and horizontal directions. Obviously, this model is multipurpose in 
real life, which has extensive application value. We refer the problem to state the general 
recursion-transform method. 

2.1 Overall train of thought and design 

We consider Ao is the origin of the rectangular coordinate system, then the left edge act 
as Y axis and the bottom edge as X axis as shown in Fig.l. Denote nodes of the network 
by coordinate (x,y). Assuming the electric current J is constant and goes from the input 
di(xi,yi) to the output d 2 {x 2 , y-i)- Denote the currents in all segments of the network as 
shown in Fig.2. The resulting currents passing through all m +1 row horizontal resistance of 
r respectively are: I^j., , • • • I^'l, (0 < k < m); the resulting currents passing through 

all n+l row vertical resistance of r'o respectively axe:lj^\ lj?\ lf \ ■ ■ ■ ij m ' ) , (0 < k < n). 

To find the resistance R m xn(dii d 2 ) we can use Ohm’s law to realize it. All nodes on 
the bottom boundary can be uniformly called O as the resistor of the bottom boundary are 
all zero. The voltages between the two nodes O and d\ , and the two nodes O and d- 2 are 
respectively 

2/1 2/2 

U m xn(0,d 1 ) = ro y £ I £ and U mxn (0, d 2 ) = r 0 £ /«. 

i= 1 i= 1 

Using Ohm’s law, we have 

R m xn(di,d 2 ) = r j (y/M - E+) (1) 

How to solve the current parameters /W and /W is the key to the problem. We are going 
to solve the problem by modeling matrix equation in terms of Kirchhoffs law. 


4 




FIG. 2: The sub-network model containing the directions of currents 


2.2 Modeling the matrix equation 


A segment of the rectangular network is shown in Fig.2. Using Kirchhoffs law ( KCL and 
KVL ) to study the resistor network, the nodes current equations and the meshes voltage 
equations can be achieved from Fig.2. We focus on the four rectangular meshes and nine 
nodes, this gives the relation, 

iffi = (2 + Miff - «ff - iffi, < = i 

iff, = (2 + 2fc)iff - hit" - hit 1 ' - iff,. 1 < * < m 

iff! = (2 + h + h 3 )lt - hit-" - 4 -!. i = m. (2) 

where hk = r^/Vo and h = v/tq. Equation (2) can be written in a matrix form after 
considering the input and output current 


Ifc+l A m xmh Ifc—1 


where I/, : and H r are respectively m x 1 column matrix, and reads 

T _ f r(l) r(2) r(3) _ _ _ r( m )]T 

- L fc L J k i 1 k i 1 k ^ i 1 k J ) 

(h^l)i T ^i,yi+l) and (-^ 2)1 ^(^*,2/2 


(3) 

(4) 

(5) 
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where [ ] T denotes matrix transposes, and (Hk)i is the elements of H r with the injection 
of current J at di(xi,yi) and the exit of current J at c/ 2 (x 2 , 1 ) 2 )■ And 


' 2 + h -h 0 ••• 0 N 

—h 2 + 2 h —h 

0 • • • — h 2 + 2 h —h 

y 0 • • • 0 —h 2 + h + h 3 j 


( 6 ) 


Next let us consider the boundary conditions of the left and right edges in the network. 
Applying Kirchoff’s laws (KCL and KVL) to a single mesh adjacent to each of the left and 
right boundaries we obtain two matrix equations to model the boundary currents, 


Ii — [A m — (2 — /?i)E]I 0 and /?,il 0 + I 2 — A m l!. (7) 

In— 1 = [A m (2 h 2 )E]I n and h 2 I n T In— 2 = A m In— i- (8) 

where E is the m x m identity matrix, matrix A m is given by (6) . 

Above Eqs.(3)~(8) are all equations we need to calculate the equivalent resistance of the 
m x n resistor network with three arbitrary boundaries. 


2.3 Approach of the matrix transformation 


How to solve the matrix equation (3) is the key to the problem. In this section we rebuild 
a new difference equation and solve Eq.(3) indirectly by means of matrix transformations. 

Firstly we consider the solution of (3) in the absence of an injected current, namely, 
J = 0. The eigenvalues U (i = 1, 2, • • •, m) of A m are the m solutions of the equation 

det |A m — t,;E| = 0. (9) 


Since A m is a special tridiagonal matrix which is Hcrmitianit can be diagonalized by a 
similarity transformation to yield 


p A IP l -1 — T 


( 10 ) 


where T m = diag{ti, t 2 , • • •, t m } is a diagonal matrix with eigenvalues ti of A m in the 
diagonal and column vectors of (P m ) _1 are eigenvectors of A m . 
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To find the exact eigenvalues t, t and the exact eigenvectors of the transformation matrix 
(Pm )' 1 is the key to the problem. Normally, it is difficult to find them. However, we can 
obtain the exact eigenvalues and the exact eigenvectors in the case of h :i = {0, h, 2 h}. These 
will be examples in the following section. 

We apply P. m on the left-hand side of (3) and write 

P?7llfc+1 P m A m I{; JP m ti X d/c X . (11) 

By (10) and (11) we define 

Pmlfc = or Ifc = P m ' X fc (12) 


where X m is a m x 1 column matrix, and reads 

y _ r y(l) y(2) y(*) y( m )]T 

■™-k lAfc i -^-k ) * " i ' ' ' i -^-k \ • 

Thus we obtain a new matrix equation after making use of (11) and (12), 


(13) 


Xfc+i — T m Xfc — Xfc_ i — JP m H x <5fc iX . (14) 

Equation (14) is much simpler than (3) as it is a linear equation easily derived. Once 
is solved that all things are easy. Next we conduct the same matrix transformation used in 
(11) to (7) and (8). Applying P m to (7) and (8) on the left-hand sides, we are led to 


A'! 1 ' = («i + /i, - 2)-C 

and + A)*- 

’ =. 

(15) 

= (U + - 2)4« 

and / 1 2 XW+XU 

) j. y (*) 

-2 ~ 

(16) 


According to the above equations we therefore obtain and X^ after some algebra and 
reduction the two solutions needed in our resistance calculation (1). 


2.4 Derivation of the resistance formula 

According to Eq.(l), the key currents and must be calculated for deducing the 
equivalent resistance R mxn (di, d 2 ). By (12) we have l k = P^X*., then we could easy obtain 
/W and /W . Finally, substituting /W and into (1), we therefore obtain the resistance 

Rmxnydli do ). 
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Obviously, the above method of computation is simple and clear, and is easy to under¬ 
stand, which has good theoretical significance and application value for the science and 
technology. As applications, we give seven new formulae by the following three examples in 
terms of our R-T method, which have not been solved before. 

3. THREE APPLICATION EXAMPLES 

In order to facilitate the expression of the equivalent resistance of resistor network, we 
define variables F^'\ /3^'j,, G$, S k: u and t t by 

jf = (A? - A*)/(A,; - A,), Aff = F®, - if, 

P& = [AF® + (ft, - lJAJ’W.JtAJ®,, + (ft 2 - 1)AF®„.J, 


G® = + (fe - l)(Ai - l)k-i + (A. + - 2)F®. (17) 

= r k /ro(k = 1, 2, 3), h = r/r 0 ; = sin((/ fc 6' i ), (18) 

U = 2(1 + h) — 2h cos 9i , (i = 1, 2.3 • • •, m), (19) 

where AjA, = 1, and A*, A,; are for later uses expressed as 

A* = 1 + h — h cos 6i + \j(1 + h — h cos 9i) 2 — 1 (20) 


Note that 0i is different in different cases, which will be determined by the different value of 
h 3 . That is to say 9 t depends on h 3 . 

3.1 Example 1 : The main result in the case of h 3 = 0 

When h 3 — 0 , from (9) and (10) we have 0 t — (i — l)7i/m, (i — 1, 2, • • •, m) ( see [23] ) , 
and 

1 l/y/2 1/V2 ••• l/y/2 

cos(l — |) 6*2 cos (2 — |) 0 2 • • • cos (m — |) 6* 2 

Pm 

^ cos(l - |)6* m cos (2 - \)6 m ■ ■ ■ cos (m - \)9 m 
The inverse matrix is (P m ) -1 = ■^■[P,, n ] T , where [ ] T denotes matrix transpose. 
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The equivalent resistance between any two nodes di(xi,yi) and d 2 (x 2 ,y 2 ) in the m x n 
resistor network can be written as (n > 1 ) 


R 


mxn 


= h 1 h 2 (y 2 -yi) 2 

m[(n - l)hih 2 + hi + h 2 ] 
r 0 yk — 2fh,2S U S 2 , + ^2,2<S , 2,i 

+ m hi (1-cos 0i)G$ 


( 22 ) 


where S kii = sin (y k 0i), 9 { = (i - l)ir/m. 

In particularly, when n —> oo,Xi,x 2 —> 00 with X\ — x 2 finite, we have result in semi¬ 
infinite case, 


R r , 


>(di, d 2 


r ™ Sf + SI - 2A^ 2 ~ ji| 5 , i5 , 2 
m i =2 \J(l + h — h cos 0 t ) 2 — 1 


(23) 


3.2 Example 2 : The main result in the case of /13 = h 


When = h , from (9) and (10) we have 9 { = (2 i — l)n/(2m + 1), (i — 1, 2, • • •, rri) ( see 
[24] ), and 

^ cos(l — \)6i cos(2 — \)9\ ■ ■ ■ 

cos(l - |)0 2 cos(2 - |)0 2 • • • 

Pm 

y cos(l - \)9 m cos(2 - 1 )9 m ■ ■ ■ 

The inverse matrix is (P m ) _1 = ^y I [P m ] r . 

The equivalent resistance between any two nodes di(xi,yi) and d 2 (x 2 , y 2 ) in the m x n 
resistor network can be written as (n > 1 ) 


cos (m 


cos (m — \) 9 \ \ 

\)9-2 

\)9 m j 


cos (m 


(24) 


^mxn(^l) ^ 2 ) 


2r„ » /31:1s, 2 , - 2/3«S M S 2lj + /3gs 2 , 


E 


kb 


(25) 


2 m + 1 i=1 (1 - cos (9j)G( 

where = sin (y k 9i), 9 t = (2 i - l)7r/(2m + 1). 

In particularly, when n —> oo,X\,x 2 —> 00 with x\ — x 2 finite, we have result in semi¬ 
infinite case, 


„ ,, ,, 2 r ” S 2 + SI - 2A1“-VS,S 2 

%xool“i,«2j — 7T5TTT 


2m + 1 i=1 ^(l + h — h cos 6^) 2 — 1 


(26) 
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3.3 Example 3 : The main result in the case of h 3 = 2 h 


When h 3 = 2 h , from (9) and (10) we have 6*, : = (2 i — l)n/(2m), (i — 1, 2, • • •, m) ( see 
[25] ), and 


( cos(l — \)9i cos (2 — \)9\ 
COs(l — |)6*2 COs (2 — 1)6*2 


cos (m 


h)»- N 


2^1 
cos(m — 5 ) 6*2 


(27) 


\ cos(l - \)9 m cos (2 - \)9 m ■ ■ • cos (m - |) 6 * m J 

The inverse matrix is (P m ) _1 = ^-[P m ] T . Note that although (24) and (27) have the same 
structure, but in which the elements 6 *; are different from each other. 

The equivalent resistance between any two nodes di(xi,yi) and < 62 (^ 2 , 1 / 2 ) in the m x n 
resistor network can be written as (n > 1 ) 


„ , . , , _ r„ A /3gSf, - 20®S u $u + 

-hmxnfOl, «2j ~ ;— 2^ - 


( 28 ) 


(l-cosSiJG? 

where S kii = sin (y k 9i), 9 { = (2 i - l) 7 r/( 2 m). 

In particularly, when n oo,X\,X 2 00 with x\ — X 2 finite, we have result in semi 
infinite case, 


„ m c2 I c2 o \ \X 2 —Xi\Q q 
r> (a a \ — r ST' "h ^2 *0102 

%x oofni, a 2) — — ; 

m i =1 y (1 + h — h cos 9 t ) — 1 


( 29 ) 


3.4 Example 4 : A general formula of the infinite resistor network 

when m,n 00 , but X\ — X 2 and yi — 1 J 2 are hnite, by (23), (26) and (29) we convert 
sums to integral by means of the limit, 

r r* 1 - A|f 2 ~' ri1 cos(j/i - 1 / 2 ) 6 * 


-6^00x00(^1, (62) 


77 J(l + h — hcos9) — 1 


(30) 


where Ae = 1 + h — h cos 9 — J (1 + h — li cos 0) — 1 . 


Please note that (30) is from different cases of h 3 = {0,h, 2h}. However, we will fully 
understand the fact since the resistance of the infinite resistor network has nothing to do 
with the boundary conditions. Thus formula (30) is a general formula of the infinite resistor 
network with four arbitrary boundaries. 
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4. APPLICATIONS OF RESISTANCE FORMULA (25) 


The above formulae obtained in this paper are too complicated to understand, we will 
give several special cases of (25) to understand its application and meaning. When h\ and 
h 2 , m or n are special number values, we have the following special cases. 

Case 1. When d\ = (xi,0) is at bottom edge, and d 2 = (x 2 , y 2 ) is an arbitrary node, from 
(25) we have 


^ 2 ) 


2 rp ™ sin 2 (y 2 6> }; ) 

2 m + 1~[ \G$ 7 1 - cos 90 


(31) 


Case 2. When d\ and d 2 are both on the left edge at d\ = (0, 2 / 1 ) and d 2 = (0, y 2 ), from (25) 
we have 


R u ,M- 2r ‘ y-/'Aei i) + (A2-l)Af’i'2 1 '|(S 1 , j -S 2 ,) 2 

n mxn\ a l, fl 2 J — o. , 1 2^ 


2 m + 1 V ) 1 - cos 0j 

Cases 3. When n —> 00 , with 04 = x 2 = 0, from (25) we have 


-Rmxoo({0, 2/l}, {0, |/ 2 }) 


2 ri 


E 


x% 1 \ {Si,i — S2,iY 


2m + l ~{\Xi + h\ — 1J 1 — cos 0, 
Cases 4. When m,n 00 , with x\ = x 2 = 0, from (33) we have 


(32) 


(33) 


■R'OOXOO ({0 ,yi},{0,y 2 }) 


n r (A e - i)[i - cos(yi - y2)0\ de 

7T Jo (\g + hi — 1)(1 — COS 9) 


(34) 


Case 5. When d\ = (0,r/i) is at the left edge, and d 2 = (n, y 2 ) is at the right edge, from (25) 
we have 


R mX n({0,yi},{n,y 2 }) 


Or m 
Zr 0 y- 

2 m + 1 “ 


hia^nSli + h 2 ai in S%j — 2hih 2 Si^S 2 ^ 

(1 — cos &i)G$ 


(35) 


where a^ n = AF® + (h k - 1)A F®i . 

Case 6 . When d\ and d 2 are both on the same horizontal axis at d\ = (xi, y) and d 2 = (x 2 , y ) 
, from (25) we have 


Rm’X.nidli d 2 ) 


2 r n 


01,1 — 201,2 + 02,2 

2m+l^V (1 — cosQYjGrJ ) 


E 


Ein 2 (y 6 >i 


(36) 


Case 7. When h 1 — h 2 — 1, the network degrades into a fan network [24], from (25) we have 


2 ftl " iSs, 2 , - 2 

2m + 1 h (l-COS e,)Fnh 


(*) 


}(0 


^mxn(^l) ^ 2 ) 


(37) 
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where /dg reduces to /dg = AFgAigg. Note that the fan network has been resolved by 
[24] in 2014, formula (37) is completely equivalent to the result in [24], Obviously the result 
in [24] is our special case in formula (25). 

Case 8. When = 0, from (25) we have 


Rr, 


i(di, d 2 ) — 


4r - /dg Si - 2 P§S lti S2,i + 


E 


(38) 


2 m + 1 ^ A i# + (hi - 1) A Eg 

where /dg is defined as /dg = [A F® + (hi - ljAFg^igg. 

The above eight cases are discovered for the first time, which shows our formulae (22), 
(25) and (28) have good application value and can deduce a lot of new results. 


5. SUMMARY AND DISCUSSION 

This paper perfect the infantile R-T method to be a complete theory, which can derive the 
general resistance between any two nodes in the nearly m x n resistor network with arbitrary 
boundaries. As applications of the R-T method, we give three examples to illuminate the 
usefulness on calculating resistance of a resistor network, and seven general formulae of the 
resistance between any two nodes in an m x n resistor network are given by the R-T method, 
which is mainly composed of the characteristic roots, and is simpler and can be easier to use 
in practice. Our three examples can’t be solved by other methods such as Greens function 
technique and the Laplacian approach. This shows the R-T method is very interesting and 
useful to solve the problem of complicated resistor network. 

In particular, our results (22), (25), (28) and (30) are four general formulae, from them 
we can deduce many simple results. Such as when we take special values of h\ — hi — 1, 
the general formula (25) immediately reduce to result in [24], Formula (30) is a new general 
result of the infinite resistor network with arbitrary boundary, which is much more concise 
and meaningful than the result in [16] and [17] since it is expressed by a single integral rather 
than a double integral. 

In addition, a very important usefulness is that the R-T method can be extended to 
impedance networks, since the Ohm’s law based on which the method is formulated is 
applicable to impedances. In addition, the grid elements r and r 0 can be either a resistor or 
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impedance, we can therefore study the m x n complex impedance networks, such as 


• r i ^ • t i R « 

r = + . , r 0 = icaL 0 + , • 

1 + icjO it 1 + icjOoJto 

Based on the plural analysis, the formula of the equivalent impedances of the m x n RLC 
network with multivariate boundary can be obtained. 
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